Radiation from a D-dimensional collision of shock waves: 
a remarkably simple fit formula 
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Recently [l| we have estimated the energy radiated in the head-on collision of two equal D- 
dimensional Aichelburg-Sexl shock waves, for even D, by solving perturbatively, to first order, the 
Einstein equations in the future of the collision. Here, we report on the solution for the odd D case. 
After finding the wave forms, we extract the estimated radiated energy for D = 5, 7, 9 and 11 and 
unveil a remarkably simple pattern, given the complexity of the framework: (for all D) the estimated 
fraction of radiated energy matches the analytic expression 1/2 — 1/-D, within the numerical error 
(less than 0.1%). Both this fit and the apparent horizon bound converge to 1/2 as D ^- oo. 

PACS numbers: 04.50.-h, 04.60.Bc, 04.30.Db 



I Introduction. A collision of two point-like particles 
at trans-Planckian centre of mass energies is thought to 
form a black hole, with an associated emission of gravi- 
tational radiation ^2, 3]. Moreover, as argued by t'Hooft 
[J|, this process should be computable by classical Gen- 
eral Relativity (GR). Understanding the physics involved 
is of conceptual interest, as a probe of the non-linear, 
non-perturbative, dynamical regime of GR, and perhaps 
even of phenomenological interest, if the current or fu- 
ture generations of particle colliders, or of cosmic ray 
detectors, can probe trans-Planckian energies. This is a 
conceivable scenario if the fundamental Planck scale is 
much lower than the traditional four dimensional one, as 
suggested in the context of D dimensional gravity (see 
[3, Sec. 4, for a recent review). 

One important property of this process is the inelastic- 
ity e of the collision: the percentage of the initial centre 
of mass energy radiated away. An upper bound for e was 
provided in [6| by an apparent horizon argument. This 
method uses no information inside the future light cone 
of the collision and can be improved by computing the 
geometry therein. A method to do so which, albeit per- 
turbative, carries information about a non-perturbative 
process, was devised by D'Eath and Payne 7-9]. This 
was recently extended to D dimensions, to first order, by 
some of us [l] . Here, we report on the case of odd D. We 
find an extra contribution, not taken into account in [l| 
(and absent for even Z?), due to the different structure 
of the odd- 1? Green's functions. Then, inspection of the 
result for 4 < Z? < 11, unveils a remarkable feature: the 
inelasticity obtained by the (technically and conceptually 
involved) D'Eath and Payne method matches, within the 
numerical precision of our method (error smaller than 
0.1%), the simple formula 
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Extrapolating this fit for large D asymptotically matches 
the apparent horizon bound. 

The inelasticity obtained for 4 < D < 11, together 



with the apparent horizon (AH) bound are shown in the 
following table (in percentage of centre of mass energy) : 
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4 


5 
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8 


9 


10 


11 


AH bound 


29.3 


33.5 


36.1 


37.9 


39.3 


40.4 


41.2 


41.9 


^Ist order 


25.0 


30.0 


33.3 


35.7 


37.5 


38.9 


40.0 


40.9 



II Wave forms in odd dimensions. The first order 
formalism of D'Eath and Payne was discussed in detail 
and generalised to higher D in 1] . Two equal Aichelburg- 
Sexl shock waves collide head on in D dimensions. The 
inelasticity of the process e can be expressed as (we refer 
to PJ for all details) 
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where the limit is selecting a radiation extraction point 
far away from the collision {r -^ oo) and along the colli- 
sion axis (0 — > 0); the wave form E,y at the space time- 
point V with null coordinates u, v and at a distance p 
from the symmetry (collision) axis is 



E^y{u,V,p) = - 
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where Av is selecting the events that support the radia- 
tion observed at V. 

The main difference between the odd D and even D 
case in Eq. (|31) is the fractional derivative of the delta 
function denoted by its exponent. The fractional deriva- 
tives of delta functions have support not only at the zeros 
of their argument but also for positive argument. This 
is related to the well known property that the retarded 
Green's function in odd dimensions has support not only 
on the past light cone, but also inside the light cone [10| . 
In other words, odd dimensional flat spacetime behaves 
like a dispersive medium. 
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FIG. 1. Wave Forms for D odd: The panels contain wave form curves for the radiation signal seen by an observer close to the 
a:xis for various (large) r as a function of time. The horizontal axis coordinate has been rescaled and shifted so that the times 
for the first and second optical rays coincide for the different curves. 



In [l[, the procedure followed for even D was to per- 
form (D — 4)/2 integrations by parts over x in Eq. ([3]), 
as to obtain a delta function and perform the x integra- 
tion completely. This procedure constrains the domain 
of integration of p' to 
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where U ^ t + 2rsm\9/2), T = t + 2r cos^{e/2) - p^ /U 
and T, r, 9 are the usual retarded time, radial and polar 
angle coordinates (see |l!|). For odd D the procedure is 
similar, except that after M = [{D — 4)/2] integrations 
by parts there is a fractional delta function of order 1/2, 
j(i/2)_ For this case, the x-integration of Eq. (B.14)of[l] 
is in fact more intricate. A careful integration by parts 
shows that the result is 
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(5) 
where now V' = V U Pextra, with 2?extra '■ Xi, < —1 , and 
the polynomial P'^*^"'"^^(a;*) factor in the domain 2?extra 
is replaced according to 
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This extra term turns out to be crucial to obtain the 
correct late time tail of the wave forms. 

A selection of wave forms is presented in Fig. [T] These 
were generated using the same numerical strategy as 



in Tr], with the extra term. We represent wave forms 
which have been rescaled by the relevant time scale for 
the problem, such that they start at retarded time ri and 
peak at T2 . Such a time scale. At = T1 — T2, is interpreted 
in the geometrical optics limit. For a (far away) observa- 
tion point not at the symmetry axis, a first ray arrives at 
the retarded time ri (corresponding to the beginning of 
the burst of radiation); then, a second ray arrives at T2, 
corresponding to the optical path that crosses a caustic 
at the axis before entering the curved region and hitting 
the observer (cf. Fig. 4 in [l|). 

Observing the curves in Fig. [T] one finds some similar- 
ities and differences with the even D case. As for the 
similarities, the number of oscillations in the wave forms 
increases with D, with one more zero for each (from left 
to right in Fig. [T|). Concerning the differences, the peak 
of radiation corresponding to the second optical ray is no 
longer singular, albeit becoming more pronounced as we 
increase r. The tails to the right of this peak are non- 
zero but integrable, since they decay as a power law. We 
have checked that the integrable tails are obtained from 
a cancellation of a non-integrable contribution from V 
and the contribution from I?oxtra- 

The inelasticity factor is also extracted numerically, 
through the double limit in Eq. ([2]). For that purpose, 
we plotted the right hand side of Eq. ^, before taking 
the limit, as a function of r in Fig. [51 for several small 
9 = TT — 9 angles. The most precise fit is extracted with 
9 = 0.01. Similarly to even D, the result is extracted 
with a relative error smaller than 10"'^. 
Ill Discussion. The method presented in [1] is tech- 
nically involved, both analytically and numerically. It is 
quite reassuring that one can obtain the results for odd 
D by the same method, fitting appropriately in the win- 
dow bracketed by the neighbouring even D values, even 
though they are obtained from integrating very differ- 
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FIG. 2. Limiting fractions: The panels contain, for each D, curves of e{r) for some values of 6 which can be used to extract the 
limiting eiat order- Some asymptotic curves which fit the numerical data to a high accuracy are indicated in each panel. The 
best estimates (relative error less than 10"'^) are indicated by the constant terms in the asymptotic fits of the red solid curves. 
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FIG. 3. Apparent horizon bound (blue dashed line and points) 
compared with the first order result (red points) and our fit - 
Eq. Il} - matching the result perfectly (red solid curve). 



ent polynomial functions. This strongly legitimates the 
method we have used. 

That the final result of this method fits, within nu- 
merical error (smaller than 0.1%), with the simple for- 



mula given by Eq. ([T|) is truly remarkable. This agree- 
ment is exhibited in Fig. [3l where the apparent horizon 
bound is also displayed. It suggests the existence of a 
simpler physical or mathematical argument to derive the 
inelasticity in this process. This certainly deserves fur- 
ther investigation and motivates the study of higher or- 
der perturbation theory for this type of processes. It is 
worth noting that, in second order perturbation theory, 
the matching conditions between the exact solution de- 
scribing the two shock waves prior to the collision and 
the perturbative method are exact. Moreover, in D = 4, 
the second order result for the inelasticity (16.3% ^) 
agrees with the value obtained in the high energy col- 
lision of two black holes in numerical relativity, within 
the numerical error (14 ± 3% [ll|). Another suggestive 
fact is the convergence of both our fit and the apparent 
horizon bound to 1/2 as I? — >■ c», but its significance, or 
if it will hold in higher order perturbation theory, is yet 
to be unveiled. 
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